An Updated Numerical Analysis of eV Seesaw with Four Generations 
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We consider the so-called "eV seesaw" scenario, with right-handed Majorana mass Ada at eV order, 
extended to four lepton generations. The fourth generation gives a heavy pseudo-Dirac neutral 
lepton, which largely decouples from other generations and is relatively stable. The framework 
naturally gives 3 active and 3 sterile neutrinos. We update a previous numerical analysis of a 3+3 
study of the LSND anomaly, taking into account the more recent results from the MiniBooNE 
experiment. In particular, we study the implications for the third mixing angle sin'^^ia, as well as 
CP violation. We find that current data do not seriously constrain more than one sterile neutrinos. 

PACS numbers: 14.60.Pq, 14.60.St 



I. INTRODUCTION 

Neutrino data hint at the possibility of something be- 
yond three massive, mostly active neutrinos. The LSND 
resuh of P{D^ -> = 0.264% [H can be explained if 
there exist one sterile neutrino, with Am^ ~ 1 eV^ with 
respect to other neutrinos. This lead to the so-called 
"eV seesaw" scenario 0], where the right-handed neu- 
trino Majorana scale Mji is taken to be of eV order. In 
2006, the eV seesaw scenario was extended to four lep- 
ton generations ^ . The fourth generation gives rise to a 
heavy pseudo-Dirac neutral lepton, which largely decou- 
ples from the other generations, and is relatively stable. 
The framework gives naturally 3 active and 3 sterile neu- 
trinos (3-1-3). 

Motivated by the interpretation of the i/g excess ob- 
served by LSND in terms of J>p — ^ Pe antineutrino oscilla- 
tions, the MiniBooNE (MB) experiment was constructed. 
In 2007, the MiniBooNE collaboration presented [l] their 
first results of a search for i^^ i/p neutrino oscillations, 
which did not support the LSND excess, but some excess 
at low energy was not fully understood. Taking this re- 
sult into account, Maltoni and Schwetz Q found that the 
3-1-2 scheme can fit both LSND i/^ data and MB 

Ve data, and at the same time account for the excess 
of low energy events in MB data. More recently, the Mini- 
BooNE collaboration reported their latest data, updating 
neutrino results @ (including the low energy region), as 
well as the first antineutrino results 0, and first results 
from the off-axis NuMI beam observed in the MiniBooNE 
detector Is]. Motivated by these new results, the authors 
of Ref. [01 re-examined the 3-1-1 and 3-(-2 global fits to 
the short-baseline (SBL) data. They found that the 3+2 
oscillation hypothesis provides only a marginally better 
description of all SBL data over the 3+1 oscillation hy- 
pothesis. In addition, a 3+1 fit to all antineutrino SBL 



data (MiniBooNE(i?), LSND, KARMEN [1(|, Bugey [11 
and CHOOZ [13]) can yield 86% x^-probability and high 
compatibility. 

In this paper we update the previous numerical anal- 
ysis of the 3+3 study of the LSND anomaly, taking into 



account MB data. In particular, we consider the implica- 
tions for the third mixing angle, and CP violation. The 
outline is as follows. In Sec. II we update the previ- 
ous numerical solutions by inputting the latest neutrino 
oscillation parameters from global data. In Sec. Ill we 
add the CP violation phases to the leptonic mixing ma- 
trix and investigate the effect of phase factors, and give 
a short summary in Sec. IV. 



II. UPDATED NUMERICAL ANALYSIS 
A. Updated Fit 

In Ref. 0, Maltoni and Schwetz discussed different 
types of oscillation frameworks that include one or more 
sterile neutrinos at the eV scale, by considering the global 
short-baseline (SBL) experiments data which included 
the MiniBooNE first result Q. The analysis built on 
an earlier study which showed that MiniBooNE (2007), 
LSND, and the null appearance experiments (KARMEN 
and NOMAD [l^) are compatible under a 3+2 ster- 
ile neutrino oscillation scenario with large CP viola- 
tion. Recently, in light of recently published results 
from the MiniBooNE experiment [a-Q, Karagiorgi et 
al. re-examined sterile neutrino oscillation models. 
They found that, with the addition of the new Mini- 
BooNE data sets, a 3+2 oscillation hypothesis provides 
only a marginally better description of all short-baseline 
data over a 3+1 oscillation hypothesis. On the other 
hand, fits to antineutrino-only data sets, including ap- 
pearance and disappearance experiments, are found sig- 
nificantly more compatible, even within a 3+1 oscilla- 
tion scenario. We shall therefore use the best-fit val- 
ues for mass splittings and mixing angles obtained from 
3+1 fits to all antineutrino SBL data in Ref. @, and the 
corresponding projected value of the LSND probability 
P(v^ Pe) = 0.2334% [13. 

We shall not recount the effective 3+3 formalism that 
arises from the eV seesaw model with 4 generations, 
which we refer to Ref. Our purpose is to assess how 



TABLE I: Best fit values, la and 3a intervals for three-flavor 
neutrino oscillation parameters taken from Ref. [isj . 



TABLE III: Parameter values of the best fit point, where the 
parameters are defined in Ref. 0]. 



BEST FIT 



Icr 



3o- 



Amii (10-5 eV^) 
Amil (10-3 eV^) 



sm 012 
sin^ 023 



sm fi3 



7.59 7.39 - 7.79 6.90 - 8.20 

+2.47 +(2.35 - 2.59) +(2.10 - 2.84) 

-2.36 -(2.29 - 2.43) -(2.00 - 2.72) 

0.32 0.30 - 0.34 0.28 - 0.37 

0.46 0.41 - 0.53 0.34 - 0.65 

0.014 0.005 - 0.025 < 0.049 



the numerics and conclusions are affected by the new 
data. To perform our numerical analysis of a 3+3 picture, 
we build the by using the three-flavor neutrino oscil- 
lation parameters Am^j and sin^ 9ij taken from Ref. [l^ , 
which are the results of the global combined analysis done 
in the framework of the AGSS09 solar fluxes jl6l and the 
modified Ga capture cross-section in Ref. [131 ■ These 
are given explicitly in Table I. We also include the pro- 
jected value of P{Df, De) = 0.2334% [Hi, and require 
Amil 1 as in the previous analysis. The seven 
inputs are listed in the second column in Table II. 

As there arc many more parameters, one cannot make 
a proper fit. so we just seek to minimize the by varying 
the model parameters. In this way, we obtain the values 
for the twelve parameters r^, Ci and Si (see Ref. 0] for 
definition) of the model at the best fit point, which are 
given in Table III. These would be further illustrated 
below. We then use these parameter values to get the 
expectation values for the seven observables, which are 
hsted in the third column in Table II. For comparison, 
the results from the 2006 study [1] is given in Table IV. 

Comparing these two tables, we can see that the differ- 
ence between best fit values and expected (input) values 
of P{Df_i -+ Pe) in Table II is smaller than the difference 
in Table IV, not just for P{i>fi — )■ Pe), but also for, e.g. 



Parameter 


ei 


£2 


£3 


Sl 


S2 


S3 


BEST FIT 


0.049 


0.018 


0.039 


0.102 


0.9999 


-0.72 


Parameter 


Til 


ri2 




r-22 


'■23 




BEST FIT 


1.263 


0.852 


1.042 


1.187 


0.742 


1.223 



sin 



623- Comparing the x^i '^'^ find 



Xmin,2006 Xmin,2010 ~ 1-83, 



(1) 



where x^i„^2006 ^nd Amin,20io a-i'e the minima for best 
fit values in 2006 and 2010. One could say that the agree- 
ment between the scenario of eV seesaw with four gener- 
ations and the experimental data has improved. 



B. In Search of Sizable sin^ 613 

The plausible value for sin^ 6*13, the target of Daya Bay 
and T2K experimemnts, is of great interest. To discuss 
this, as in Ref. [1], we restrict the to the four inputs 
of sin^ Oij and P{Df^ v^)- With and given as in 
Table III, we first fix si = 0.1, the best fit value, and 
perform a fit vs S2 and S3, then iterate with fixing 
S2 = 0.9999 (s3 = —0.72) and minimize vs si and S3 
(si and S2). We find for both cases of fixing si or S3 to 
the best fit values, sin^ ^23 is quite strongly dependent 
on S2, and the value around 0.9999 is preferred. We thus 
illustrate with S2 held fixed to this value. 

In Fig. 1 we show the x^ contour plot vs si, S3. The 
three different shaded regions should not be interpreted 
as the 1(7, 2(7 and 3(7 regions, since we have fixed the 
rest of the parameters to the best fit values. But they do 
give an indication of variations around the best fit region 
under the above assumptions. The gross features are not 
too different from the 2006 study. For the lower right 



TABLE II: Seven inputs taken from Ref. [1, [3, ^ and the 

best fit values from the minimization process. 



TABLE IV: Inputs [I, [Ti| and the best fit values from in 
2006 [ll, to be compared with Table II. 





INPUT 


2010 




INPUT 


2006 




0.2334 


0.2090 




0.264 


0.15 


Amil (lO-'^ eV^) 


7.59 


7.59 


Amil (lO-'^ eV^) 


8.1 


8.1 


Amil (10-3 eV^) 


2.5 


2.5 


Amil (10-3 eV^) 


2.2 


2.3 


Amil (eV^) 


1 


1 


Amil (eV") 


1 


1 


sin^ ^12 


0.32 


0.32 


sin^ 9i2 


0.30 


0.30 


sin^ 623 


0.46 


0.46 


s\v? O23 


0.50 


0.52 


sin^ 613 


0.014 


0.014 


sin" 613 


0.0 


0.0018 
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Fig. 1: Contour-plot of vs the mixing angles si and S3, 
with Ei and Vij as in Table IV, and S2 — 0.9999 held fixed. 
The regions in different shades are only indicative, and should 
not be interpreted as the la, 2a and 3a regions, as the rest 
of the parameters are fixed at the best fit values. 



solution, si has moved closer to zero, and the strength of 
S3 has also weakened a little. For the upper left solution, 
S3 is little changed, but si has strengthened considerably, 
moving the solution more into the corner. 

We plot in Fig. 2 the four quantities PiPfj. S'e), 
sin^6'i2, sin^ 6*23 and sin^ ^13 vs si and S3, for the solu- 
tion on the lower right of Fig. 1. The same is plotted in 
Fig. 3 for the upper left solution. Again, and are 
fixed as in Table III, and S2 is held fixed at 0.9999. The 
dependence on S1-S3, or shape, for /'(P^^ i>e) differs 



P(v^ ^ Ve)% sin^(ei2) 




Fig. 2: P(i'n — !> i'e), sin^6'i2, sin^ 623 and sin^ 6^13 vs si and 
S3, corresponding to the lower right solution in Fig. 1, with 
ti and Tij fixed as in Table III, and S2 fixed at 0.9999. 



somewhat from the 2006 ease, with the lower right solu- 
tion now allowing larger values. We see that P{v^ Ve) 
can reach 0.2% and sin^ 612 is well within range, and 
likewise for sin^6'23. However, to push sin^ 6'i3 beyond 
0.025, sin^ 623 would start to wander away from maximal 
mixing of 0.5, and values at ^ 0.4 or 0.6 has to be toler- 
ated. The situation is not much changed from 2006, but 
we note that the slightly more promising situation for 
sin^ ^13 can be traced to the more positive input value in 
Table II. which now docs include SBL data. 



C. Neutrino Mass Hierarchy 

The sign of Amgj^ is undetermined by present data 
(e.g. MINOS experiment). That is to say, the neutrino 
mass hierarchy is not known at this time. We use eV see- 
saw with four generations scenario to see if we can gain 
any access to the neutrino mass hierarchy. For normal 
mass hierarchy, the sign of Amgj^ is positive and the best 
fit value is 2.47 x lO"'^ eV^, and the numerical analysis 
has been performed in the previous section. For inverted 
mass hierarchy, the sign of Am|]^ is negative with best fit 
value —2.36 x 10"'^ eV^. To perform the numerical anal- 
ysis for inverted mass hierarchy, we use the best fit values 
given as in Table III but only replacing that of Amgj^ by 
—2.36 X 10""^ eV^ to work out minimum. Comparing 
with the minimum for normal mass hierarchy, we find 

Xmin, inverted Xmin, normal 1.(4, (2) 

in obvious notation. It seems that the agreement of eV 
seesaw with four generations scenario with normal mass 



P(v,, ^ Ve)% sin^(ei2) 




Fig. 3: Same as Fig. 2, but corresponding to the upper left 
solution in Fig. 1. 
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hierarchy is sUghtly better than the agreement with in- 
verted mass hierarchy. We caution that our study is not 
a true fit, as there are too many model parameters. 

III. CP VIOLATION 

So far, we have not considered CP violation effect, and 
all phases in the leptonic mixing matrix were set to zero 
for simplicity, just as in the previous paper Q . As de- 
fined in Ref. Q, U' and U" are the rotation matrices 
which diagonalize the neutrino and charged-lepton mass 
matrices, respectively, and the full leptonic mixing ma- 
trix is [/ = U"U'. Although introducing CP violation 
phase to the leptonic mixing matrix adds further to the 
already many parameters, there is some motivation to 
study CP violation effect, in part because there are both 
neutrino and antineutrino oscillation data, and they ap- 
pear to be somewhat different. We shall therefore give 
two examples for illustration. 

A. Refit with CP Violation Pliase 

As charged Icpton masses are generated by the Higgs 
mechanism, similar to quark masses, we take the left sec- 
tor rotation matrix U" to be the usual form 



/ 


C1C3 


S1C3 


S3e^'^ 












— S1C2 — CiS2S3e'* 


C1C2 - SiS2S3e"' 


S2C3 













S1S2 - CiC2S3e'* 


-C1S2 - SlC2S3e"' 


C2C3 






















1 






















1 





V 

















1/ 



(3) 

where 6 is the CP violation phase and placed in the 13 
element of U" . Because of the presence of S, the formu- 
las for De) would no longer be equal to that of 
Pivf^ Ve). We use this new U" matrix with CP viola- 
tion phase (5, and the best fit values given as in Table II 
for inputs, and redo our numerical analysis. After mini- 
mization, we get the new minimum, i.e. Xmin cp ■ 
find that introducing CP phase 5 leads to the relative 
improvement of the fit of 

Xmin,2010 - Xmin,CP = 0-004, (4) 

which is negligible. The quality of the global fit does not 
improve at all by introducing CP violation phase. But 
this is somewhat trivial, since Table II tacitly assumes 
CP conservation. 



B. An Extreme Example of CP Violation 

The updated MiniBooNE result indicates an 

excess of events at low energy, but no oscilla- 
tions at the L/E ^ 1 m/MeV of the LSND excess, which 
is for antineutrinos. Furthermore, neither NUMI Q nor 



TABLE V: Same as Table II, but adding a null — !> Ve 
oscillation probability as input, together with the expectation 
values from the minimization process. 



INPUT 





0.2334 


0.1423 







0.048 


l\ml^ (10-5 eV^) 


7.59 


7.59 


Amii (10-3 eV^) 


2.5 


2.5 


Am|i (eV^) 


1 


1 


sin^ 6*12 


0.32 


0.32 


sin^ 6*23 


0.46 


0.43 


sin^ 6*13 


0.014 


0.017 



NOMAD found evidence for — > Ve oscillations. 
Thus, we consider P(i/^ v^) = 0% at L/E = 1 m/MeV 
as a new additional input, to illustrate an extreme sit- 
uation where there is a striking difference between the 
behavior of i'^ — !■ Ve and Dfj, — >■ De oscillations. We 
stress that this is only an illustration, since neither the 
LSND effect for Df^ — > i>e oscillations, nor the absence of 
j/^ — > j^e, are firmly established. Furthermore, we clearly 
have too many parameters, hence we are just performing 
illustrations, not a fit. 

The new ad hoc input, together with the 7 inputs of 
Table II, arc now listed in Table V, with which we build 
the new in the framework with CP violation. After 
minimizing the x^ j we use the thirteen parameter values 
of the best fit point to get the the expectation values for 
the eight observables, which are listed in the third col- 
umn in Table V. There is now a slight tension in sin^ 623 
(though allowed), while neither the antineutrino, nor the 
neutrino oscillation probabilities can reach the input val- 
ues, but this may be reasonable. The best fit value of 
CP violation phase S is 0.377r, or of order 66°, which is 
of course rather large. Note also that the sin^ 613 value 
is now larger than in Table II, which is again reasonable, 
since a finite sin^ ^13 is a prerequisite for CP violation. 
This illustrates the efficacy of the Daya Bay and T2K 
experiments, as well as a future CPV program. 



IV. SUMMARY 

We have updated the numerics of a 3-1-3 active/sterile 
neutrino picture, which arises naturally from a four lep- 
ton generation model with eV seesaw, i.e. with right- 
handed Majorana mass scale Mn at eV order. The ma- 
jor new data are the MiniBooNE results of the past few 
years, which do not confirm the LSND D^^ — >■ result. 
Since the 3-1-3 model has too many parameters, while 
data is more or less consistent with a 3-1-1 picture, we 
take global fit results using the latter framework as input. 



4 



and use a ^ minimization to find the best values for neu- 
trino mass differences, mixing angles, and — > prob- 
ability. Not unexpectedly, the consistency of the 3+3 
picture with current data improves over the consistency 
with data in 2006. We continue to find it difficult for 
sin^ 6113 to be larger than 2.5%, while finding a slightly 
better for normal mass hierarchy over the inverted 
mass hierarchy. These, however, could be artefacts of 
not being able to do a real fit. 

If one continues to take the LSND v^^ Vf. result se- 
riously (which provides the eV scale), together with the 
consensus of null results for v^^ — > Vf, oscillations, we have 
put in the CP violation phase in the charged lepton mix- 
ing matrix, and illustrated the possibility of CP viola- 
tion. The world of large CP violation effect in — ?> Ve 
vs Ue oscillations, with large CPV phase and finite 



sin^ 013, could be a striking one. 

Note Added. Our Table I uses version 2 of Rcf. [isj . 
which is the published one. Subsequently, these authors 
have posted version 3 on the arXiv, which includes new 
results on appearance from MINOS, which leads to 
a smaller sin^ 6'i3 (best fit at 0.008). However, we have 
discussed how things will be for a lower sin^0i3, so we 
leave the paper as is. 

Acknowledgement. We thank T. Schwetz and 
G. Karagiorgi for useful communications. The work of 
WSH is supported in part by NSC97-2112-M-002-004- 
MY3 and NTU-98R0066. The work of FFL is supported 
by NSC98-2811-M-002-102. 



[1] A. Aguilar et al. (LSND Collaboration), Phys. Rev. D 

64, 112007 (2001). 
[2] A. de Gouvea, Phys. Rev. D 72, 033005 (2005). 
[3] W.-S. Hou and A. Soddu, Phys. Lett. B 638, 229 (2006). 
[4] A. Aguilar- Arevalo et al. (MiniBooNE Collaboration), 

Phys. Rev. Lett. 98, 231801 (2007). 
[5] M. Maltoni and T. Schwetz, Phys. Rev. D 76, 093005 

(2007). 

[6] A. Aguilar- Arevalo et al. (MiniBooNE Collaboration), 

Phys. Rev. Lett. 102, 101802 (2009). 
[7] A. Aguilar- Arevalo et al. (MiniBooNE Collaboration), 

Phys. Rev. Lett. 103, 111801 (2009). 
[8] P. Adamson et al, Phys. Rev. Lett. 102, 211801 (2009). 
[9] G. Karagiorgi, Z. Djurcic, J.M. Conrad, M.H. Shaevitz 

and M. Sorel, Phys. Rev. D 80, 073001 (2009); Erratum- 

ihid. D 81, 039902(E) (2010). 
[10] B. Armbruster et al. (KARMEN Collaboration), Phys. 



Rev. D 65, 112001 (2002). 
[11] B. Achkar et al, Nucl. Phys. B 434, 503 (1995). 
[12] M. Apollonio et al (CHOOZ Collaboration), Eur. Phys. 

J. C 27, 331 (2003). 
[13] P. Astier et al (NOMAD Collaboration), Phys. Lett. 

B570, 19 (2003); D. Gibin, Nucl. Phys. Proc. Suppl. 66, 

366 (1998). 

[14] We are grateful to Thomas Schwetz for providing us this 
value. 

[15] M.C. Gonzalez-Garcia, M. Maltoni and J. Salvado, JHEP 
1004, 056 (2010) |arXiv:1 001.4524v2 [hep-ph]]. 

[16] A. Serenelli, S. Basu, J. W. Ferguson and M. Asplund, 
arXiv: 0909.2668 [astro-ph.SR]. 

[17] J.N. Abdurashitov et al. (SAGE Collaboration), Phys. 
Rev. C 80, 015807 (2009). 

[18] M. Maltoni et al, New J. Phys. 6, 122 (2004). 



5 



